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SECTION  1 


INTRODUCTION 

A primary  consideration  in  the  design  of  structural  sandwich  com- 
ponents is  the  overall  stability  of  a panel  which  is  subjected  to  inplane 
loadings.  These  applied  loads  commonly  consist  of  edgewise  compression, 
edgewise  shear  and  edgewise  bending  moments,  as  shown  in  Figure  1. 

A large  amount  of  literature  is  available  concerning  the  buckling  of 

flat,  rectangular  sandwich  panels  loaded  by  uniaxial  edgewise  compression 

7 

Kuenzi,  et  al.  have  studied  the  problem  of  shear  acting  alone.  Other  inves- 
tigators have  considered  the  effects  of  compression  combined  with  shear  by 

8 9 

the  use  of  interaction  formulas  ’ designed  to  account  for  interference 
between  the  various  buckle  patterns.  The  case  of  combined  edgewise  bending 
and  compression  has  been  investigated  by  Kimel  with  good  results.  The 
results  of  a similar  study  which  considered  combined  bending,  compression 
and  shear^  exhibited  poor  convergence  behavior,  particularly  for  cases  in 
which  the  shear  loads  predominated. 

This  report  presents  a method  for  the  prediction  of  the  general 
instability  of  flat  sandwich  panels  loaded  by  combined  edgewise  bending, 
biaxial  compression  and  edgewise  shear.  It  is  anticipated  that  the  method 
discussed  herein  will  prove  useful  as  a tool  in  both  the  design  and  the 
analysis  of  structural  sandwich  components. 

1.  1 PROBLEM  DEFINITION  AND  SCOPE 

The  present  analysis  considers  incipient  buckling  of  flat,  rectangular 
sandwich  panels  subjected  to  arbitrary  combinations  of  biaxial  edgewise 


Numerical  superscripts  indicate  references  listed  in  Section  6. 
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Figure  1.  Applied  Inplane  Loads. 
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compression,  edgewise  bending  and  edgewise  shear  loadings.  Only  simply 
supported  panels  are  considered,  for  which  the  component  layers  are 
restricted  to  linear  elastic  behavior  and  small  displacements. 

The  core  is  considered  to  be  completely  rigid  in  the  normal  direction. 
The  resistance  of  the  core  to  extension  and  bending  are  assumed  to  be 
negligible;  therefore  the  strain  energy  may  be  taken  to  consist  only  of  con- 
tributions due  to  transverse  shearing  deformation.  The  transverse  shear 
moduli  of  the  core  may  be  unequal,  as  is  typical  of  many  types  of  core 
construction. 

Facing  sheets  are  assumed  to  be  isotropic  and  may  have  different 

$ 

properties  and  thicknesses.  Since  the  flexural  rigidity  o'  the  faces  is  not 
neglected,  critical  combined  loads  of  classical  plates  may  be  obtained  as  a 
special  case  (subject  to  the  Love-Kirchhoff  hypothesis). 

1.  2 METHOD  OF  SOLUTION 

In  general,  linear  analyses  of  sandwich  deformations  fall  into  two 
1 2 

categories:  formulations  in  terms  of  face  displacements,  and  analyses 

in  terms  of  functions  related  to  core  displacements.  The  two  methods  are 
theoretically  equivalent  but  permit  slightly  different  numerical  treatments 
of  the  boundary  conditions. 

The  results  presented  herein  are  based  upon  the  so-called  "tilting 
1 3 

method"  , by  which  the  displacements  in  each  layer  of  the  sandwich  are 
related  to  displacements  in  the  core.  A solution  is  obtained  by  using  a 
discretization  based  on  the  Ritz  method  to  minimize  the  total  potential 
energy  of  the  panel.  As  a means  of  uncoupling  the  flexural  and  extensional 
strain  energies,  it  is  postulated  that  a "neutral  plane"  exists  in  which  no 
inplane  displacements  occur  during  the  buckling  deflection.  The  neutral 

*The  accurate  analysis  of  severely  unbalanced  sandwich  may  require  the 
use  of  a nonlinear  technique.  See  Paragraph  1.  2. 
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r 
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plane  location  which  corresponds  to  a given  buckling  mode  is  then  treated 
as  an  additional  degree  of  freedom  in  the  process  of  minimizing  the  energy 
expression. 

It  should  be  noted  that  such  a linear  analysis  is  valid  only  under 
certain  ideal  circumstances.  The  most  important  of  these  are  the  require- 
ments that  the  structural  components  under  consideration  experience  only 
infinitesimal  changes  in  geometry  prior  to  buckling,  and  that  sufficient 
symmetry  is  present  that  any  assumptions  aimed  at  uncoupling  modes  of 
deformation  (in  this  case,  the  existence  of  neutral  planes)  are  justified. 

In  particular,  sandwich  panels  of  severely  unbalanced  face  construction 
are  likely  to  require  a more  rigorous  (i.  e.  , nonlinear)  treatment  than  is 
presented  here,  since  the  range  of  validity  of  the  assumptions  concerning 
the  neutral  planes  may  be  exceeded. 
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SECTION  2 

THEORETICAL  APPROACH 

A displacement  formulation  is  used  to  construct  the  total  potential 
energy  of  the  panel.  The  theory  upon  which  this  analysis  is  based  is  pre- 
sented in  Reference  12.  The  essential  equations  are  cited  here  for 
completeness. 

2.  1 GEOMETRY 

The  panel  configuration,  reference  axes  and  loading  conventions  are 
indicated  in  Figure  2.  The  x-y  plane  is  chosen  to  coincide  with  the  panel 
midplane.  Positive  senses  of  the  stress  resultants  and  displacements  are 
as  shown  in  the  figure. 

Figure  3 shows  the  assumed  geometry  of  deformation  in  the  xz 
coordinate  plane.  The  displacement  un  represents  an  inplane  deflection 
due  to  stretching  only,  and  is  not  of  interest  here,  since  only  the  state  of 
incipient  buckling  is  to  be  considered.  As  the  panel  buckles,  a transverse 
deflection  occurs,  accompanied  by  transverse  shearing  deformations  within 
the  core.  During  this  deflection,  it  is  assumed  that  in  some  plane  located 
at  z = e there  is  no  further  inplane  displacement,  and  that  the  rotations  of 
the  panel  in  the  xz  plane  occur  about  points  in  this  plane.  Similarly  in  the 
yz  plane,  rotations  are  assumed  to  take  place  about  the  plane  z = e^.  With 
this  assumption,  each  of  the  inplane  displacements  of  the  panel  are  expressed 
in  terms  of  the  transverse  deflection  w and  the  core  rotations  4>  and 

2.  2 BASIC  EQUATIONS 

The  assumed  geometry  (Figure  3)  permits  the  displacements  of  the 

1 2 

core  and  faces  to  be  stated  as 
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where 


and 


“fy  = V ('1)f  ey  5 1 = l>  2 • (2.  2) 


It  is  important  to  notice  that  the  neutral  plane  locations  e , e are  unknown 

x y 

parameters  to  be  determined  in  the  process  of  minimizing  the  potential 
energy. 


The  strains  are  then  found  to  be 
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and  the  stresses  by  Hooke's  law 


■j  c - z ( * , + V t ,) 

x*  i v 2 x£  £ yf' 
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r = G y 
cxz  cxz  cxz 


r - G y 
cyz  cyz  cyz 


It  is  evident  that  each  of  the  components  of  displacement,  strain  and  stress 

due  to  bending  of  the  sandwich  is  a function  only  of  the  core  rotations  </> 

and  <!>  and  the  transverse  displacement  w.  The  responses  due  to  inplane 

stretching  are  functions  of  the  displacements  u and  v . 

n n 

Notice  in  Equation  2.4  that  the  inplane  stress  resultants  applied  to 

the  panel  up  to  the  point  of  incipient  buckling  are  those  related  to  the 

stretching  deformations  u and  v ; that  is, 

n n 


N = 
x 


d/2 

I - 


-d/2 


,1-  v 


- (u  + v v ) dz 
2 n,x  n,y 


/ 2ufn  <°,.y  f v„.x,ds 


•d/2 


(2.5) 


E.i  ’ ° : 

; , j (v  + V u ) 

• i>v2  n>y  ' : n’x 


Since  the  core  does  not  resist  inplane  stretching,  these  integrals  receive 
nonzero  contributions  only  from  the  face  sheets.  Accordingly,  each  of 
the  resultants  defined  by  Equation  2.  5 may  be  thought  of  as  the  sum  of  the 
stress  resultants  in  the  two  facings;  that  is. 


N = N + N , 
x xl  x2 


N = N + 

y yi 
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y2 


N 

xy 


+ N 


xy2  * 


(2.6) 


2.  3 POTENTIAL  ENERGY 


IT  = - 

P 2 


The  total  potential  energy  of  the  panel  may  be  expressed  as 
2 

2 £ Iff  kf'xf*  rxyf*xyf  + °yf  V 
f=l  V,  L 


14 


1-2  - — 2 * 
+ — (N  w,  + 2 N , w,  w + N w,  ) 
tf  v xf  x xyf  x y yf  y ' 


dxdydz 


( r y + T y ) dxdydz 
cxz  cxz  cyz  cyz 


(2.7) 


where  f denotes  the  respective  face  sheets,  and  and  represent  the 

volumes  of  the  facings  and  the  core,  respectively.  In  the  statement  of 

Equation  2.  7,  it  is  tacitly  assumed  that  displacements  experienced  at  points 

within  the  panel  prior  to  buckling  are  negligible  (that  is,  ur  = vn  = 0).  The 

stress  resultants  N „ N ,,  N , have,  therefore,  assumed  the  character  of 
xf  xyf  yf 

initial  stresses  present  throughout  the  face  sheets  immediately  preceding 
the  time  of  buckling.  Furthermore,  the  stresses  and  the  strains  are 
to  be  interpreted  as  departures  from  the  prebuckled  equilibrium  state 


10 


defined  by  the  displacements  u and  v . It  should  also  be  noted  that  the 

n n 

squares  of  the  derivatives  of  the  inplane  displacements  have  been  neglected 
in  comparison  with  the  squares  of  the  slopes;  that  is, 


2 2 2 2 2 2 
(u  ; u.  ; v,  ; v ) <<(w,  ; w,  ) 

f,x  f,y  f,x  f,y  x y 


(2.8) 


The  displacement  formulation  corresponding  to  Equation  2.  7 is 
obtained  with  the  use  of  Equations  2.  3 and  2.4.  It  is  then  assumed,  for  the 
linearized  problem,  that  cubic  and  quartic  terms  in  the  strain  energy  may 
be  neglected.  After  making  this  substitution  and  performing  the  integrations 
through  the  thickness  of  the  panel,  the  following  potential  energy  functional 
is  obtained: 


2 E b a / 


f = 1 o o 


. 2 t / 2 l-vf  2 

'f*  f U, 

' X Z v 


+ «.  t *•  + +>  i 

fx  fy  I f x y f y xl 


2 . 2 f . 2 \ 

+ «fy'f  V ’ y +—  *•.  ) 


I t k (w,  + V W,  ) + (1-V,)<£,  W,  1 

fx  f L x xx  f yy  f y xyj 


t,  k <W*  + V r W,  ) + (1-V  , W,  1 

r f l y yy  f xx  fx  xyj 


1 3 r 2 

T t w,  + 2v  w,  w,  f w, 

3 f L xx  f xx  yy  yy 


+ i'1-V('*’xy  ]} 


dxdy 


This  assumption  is  reflected  in  Equation  2.  3. 

2 

The  definition  A^=l-v  ^ has  been  introduced  here. 
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b a 


\ f f\c  t ( + w,  + G t ( «£  + w,  ) |dxdy 

2 J J L cxz  c x cyz  c y J 7 


o o 


b a 


-7  f / [ N w/  + 2N  w,  w,  + N w, 2 1 dxdy 
2 J J L x x xyxy  yyj  7 


(2.9) 


o o 


It  is  important  to  realize  that,  although  the  resultants  N , N N are 

x xy,  y 

required  to  satisfy  equilibrium  before  buckling,  they  need  not  be  constants 
throughout  the  panel.  Hence,  the  quantities  and  N may  be  taken  to 
include  both  uniform  and  linearly  varying  compressive  forces;  that  is, 


N = N + N (y  / b ) 
x xo  xb 

C. 


^y  = N + N (x/a) 
yo  yb 


(2.  10) 


These  forms  for  the  applied  loads  satisfy  the  equations  of  equilibrium. 


N + N 

= 0 

x,x  xy,  y 

N + N 

= 0 , 

xy,x  y , y 

and  permit  the  consideration  of  edgewise  bending  forces.  The  special  case 

of  a pure  edgewise  bending  moment  (see  Figure  1)  will  be  indicated  by  the 

symbols  N „ or  N , such  that 
xB  yB 


NxsNxB(1-^/b) 


N = N u ( 1 - 2x/a ). 

y yH 


(2.  11) 


That  is,  the  linearly  varying  component  of  the  pure  edgewise  bending  load 
is  twice  the  magnitude  of  the  uniform  component,  and  opposite  in  sign. 
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SECTION  3 


NUMERICAL  SOLUTION  FOR  BUCKLING  LOADS 

The  principle  of  minimum  potential  energy  can  now  be  applied  to 
produce  the  conditions  of  equilibrium.  The  potential  energy  is 
expressed  in  terms  of  a complete  and  geometrically  admissible  set  of 
functions.  Application  of  the  principle  of  minimum  potential  energy  pro- 
duces the  approximate  minimizing  relations,  in  the  form  of  an  eigenvalue 
problem.  By  considering  successively  larger  solution  spaces,  each  of 
which  contains  the  previous  one,  monotonic  convergence  to  the  solution 
is  guaranteed. 

3.  1 DISCRETIZATION  OF  THE  POTENTIAL  ENERGY  BY  THE 
RITZ  METHOD 

The  total  potential  energy  associated  with  the  incipient  buckling  of  the 
panel  is  given  by  Equation  2.9  in  terms  of  w,  <f>  and  Admissible  assumed- 
mode functions  must  satisfy  conditions  of  continuity  and  differentiability,  as 

12 

well  as  the  imposed  boundary  conditions  for  simple  supports, 
w(o,y)  = w(a,  y ) = w(x,o)  = w(x,b)  = 0; 
therefore  an  appropriate  set  of  functions  is 


i j e 

w(x,y)=  Z Z W sin  mTrX  sin 

’ ' mn  a b 

m=l  n = l 


l l 

0(x,y)  = Z Z ♦ cos  r^TrX-  sin 

' mn  a b 


(3.  1) 


m=l  n = l 


l l 


*<x.y)=  Z £ *mn8in^C08  T* 


m=  1 n = l 


-'“•4U 


Substitution  of  the  assumed  forms  of  the  displacements  into  the  potential 
energy  (Equation  2.9)  yields  a functional  u^  which  is  expressed  entirely  in 
terms  of  the  trigonometric  functions  of  Equation  3.  1.  The  integrations 
over  the  area  may  then  be  carried  out  directly,  according  to  the  following 
formulas: 


a 


/ 


mirx  rnx  , a 

sin  sin dx  = — 

a a 2 


8 

mr 


/ 


mirx  r-rrx  , a . 

cos  cos  dx  = — 8 

a a 2 mr 


/ 


. mux  rux  , 2am  A 

sin c os  dx  = A 

a a ,2  2.  mr 

ufm  - r ) 


/ 


mux  . rux  , a 
xsm  sin dx  = — 8 


. 2 

4a  mr 


A , 


4 mr  2,  2 2.2  mr 

u (m  - r ) 


where  8 is  the  Kronecker  delta,  and  A is  defined  such  that  for  any 
mr  mr  r 

arithmetic  expression  denoted  by  X, 


!0  ; m+  r even 
X ; m+  r odd 

The  expression  obtained  for  the  potential  energy  upon  substituting  the 
assumed  modes  (Equations  3.  1)  and  performing  the  indicated  integrations 
is  given  by 

♦ The  expression  may  be  indefinite  or  definite. 
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It  is  important  to  recognize  that  since  the  parameters  «„  , a contain  the 

£x  fy 

positions  of  the  neutral  planes  (Equation  2.  Z),  they  too  are  undetermined 
parameters  dependent  upon  the  particular  buckling  mode;  that  is, 

« , = ( «,  ) 
fx  fx  mn 

fy  fy  mn 

Since  the  energy  expression  contains  cubic  terms  (such  as  a <l>2  ) 

fxmn  mn 

which  are  not  quadratic  forms  in  the  unknowns  4>  , * , W , e 

mn  mn  mn  xmn 

and  e , the  minimization  of  the  energy  in  its  present  form  will  involve 
ymn  r 

the  solution  of  a system  of  nonlinear  simultaneous  equations.  In  order  to 
circumvent  this  difficulty,  it  is  advantageous  to  define  the  additional 
parameters 

h = e <t> 

mn  xmn  mn 

k = e 4 

mn  ymn  mn 

It  is  also  convenient  to  define  a (4  x 4)  matrix  A,  a (4  x 1)  vector  B,  and  a 
constant  C corresponding  to  each  mode,  as  follows: 
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With  these  conventions  and  the  use  of  Equation  2.2,  the  potential  energy 
(Equation  3.  2)  can  be  written  in  the  form 
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(3.3) 


The  discretized  potential  energy  (Equation  3.  3)  is  a quadratic  function  of  the 

undetermined  parameters  4>  , ♦ , h , k , and  W . The  use  of  the 

mn  mn  mn  mn  mn 

principle  of  minimum  potential  energy  then  leads  to  a linear  system  of 
equations  for  the  undetermined  parameters. 


3.2  APPLICATION  OF  THE  PRINCIPLE  OF  MINIMUM  POTENTIAL 
ENERGY 

The  variational  form  of  the  principle  of  minimum  potential  energy  is 


4 it  =0, 

P 

where  rr  is  given  in  Equation  2.9.  After  discretization  using  the  Ritz 
P 

method,  the  minimum  potential  energy  principle  requires 

gp  ^p  ^1*  ^2*  * * * * Pn^  = ^ 2,  . . . . , n. 


where  the  p.  are  unknown  discrete  parameters.  For  the  present  case, 
from  Ecuation  3.  3, 
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The  unknown  parameters  may  be  normalized  with  respect  to  by  taking 
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Upon  substituting  these  definitions  in  Equations  3.4  and  3.5,  the  following 
final  system  of  equations  is  obtained: 
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(3.7) 


3.  3 ASSEMBLY  OF  EQUATIONS  FOR  SOLUTION 

It  is  clear  from  Equations  3.  6 and  3.  7 that  only  the  parameters  W^ 
must  be  considered  explicitly  in  the  solution  process,  since  Equation  3.  6 can 
be  inverted  to  obtain  the  X...  The  quantities  are  collected  in  a vector, 


MT  = l»1i-*i2> 


.,wu.w21.w22,...w 


2 i wf  r Wi2'  • • 'wU—t ■ 
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< 


The  prescribed  loads  are  assigned  constant  relative  magnitudes,  and  are 
related  to  a single  factor  which  is  sufficient  to  indicate  their  intensity; 
that  is, 


( N , N , , N , N , N ) = X ( n , n , n , n , n ). 
xo  xb  xy  yb  yo  xo  xb  xy  yb  yo 


Equations  3.  7 can  then  be  arranged  in  the  form 


[k]  {w}  + »[g] 


(3.8) 


2 2 2 
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2 2 2 
l x l lx  1 


Here  the  matrix  K is  diagonal,  with  the  i^  diagonal  entry  defined  by 
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k..  = C b.  X. 

li  mn  / j jmn  jmn 


where 


j = l 

i = X(m-  1 ) + n 


and  the  X are  solutions  of  Equation  3.  6 for  the  appropriate  m and  n.  The 
th  J 

i,  j entry  of  the  matrix  G is  related  to  the  applied  loads  terms  in 


Equations  3.  7,  such  that 
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whe  re 


i = X(m-  1 ) + n 


and 


j = i(r-l)  + s . 

Finally,  the  critical  loads  and  corresponding  mode  shapes  of  the  panel  are 
obtained  by  solving  the  eigenvalue  problem  represented  by  Equation  3.8. 
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SECTION  4 


SUMMARY  OF  RESULTS 

The  procedure  outlined  in  Section  3 has  been  implemented  in  a com- 
puter program.  Numerical  results  have  been  obtained  for  numerous  examples 
of  individual  critical  loads  from  other  sources  in  order  to  verify  its 
accuracy.  Critical  loads  have  also  been  calculated  for  various  cases  of 
combined  loadings,  and  compared  with  existing  interaction  formulas. 

4.  1 VERIFICATION  EXAMPLES 

The  results  of  typical  calculations  performed  in  the  process  of 
verifying  the  computer  program  are  presented  in  Table  1.  The  example 
cited  from  Reference  15  is  based  upon  a nonlinear,  post  buckling  finite 
element  analysis.  The  remainder  are  taken  from  linear  buckling  analyses 
similar  to  that  presented  here.  Note  that  Examples  4,  5,  6,  and  11  are 
specializations  of  the  theory  to  the  case  of  a classical  isotropic  plate. 

Agreement  with  the  previously  published  results  was  good  in  all 
cases.  Critical  loads  calculated  for  compression  agree  closely  with  those 
cited,  since  the  various  modes  uncouple,  allowing  the  equations  to  be 
solved  in  closed  form. 

* 

Pure  edgewise  bending  results  are  also  in  good  agreement  with 
the  previous  values.  Convergence  for  these  calculations  generally  requires 
no  more  than  6 terms  in  the  displacement  series  (Equation  3.  1).  The 
results  cited  from  Reference  10  apply  for  face  sheets  having  only  membrane 
stiffness.  Consequently  the  present  analysis  yields  a slightly  higher  critical 
load  for  Example  9. 


Pure  edgewise  bending  values  are  stated  in  the  form  of  Equations  2.  11; 

that  is,  the  value  of  N defines  N = N n (l-2y/b),  and  so  on. 

x Jd  X XD 
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TABLE  1:  VERIFICATION  RESULTS  FOR  INDIVIDUAL 

CRITICAL  LOADS 


No. 

Construction 

Loading 

E, 

Ez 

'i 

'2 

*1 

‘2 

G 

cxz 

G 

cyz 

a 

H 

Terms 

in 

Series 

Critical 

Loads 

— 

Ref. 

Reference 

Result 

■ 

9.  5xl06 

. 30 

.021 

19000. 

. 181 

R5B 

Nx«  -309.  39 

15 

-307.  5 

H 

Sandwich 

Compr. 

9. 5x10 

. 30 

.021 

19000. 

— 

RVHKj 

. 30 

.025 

50000. 

.25 

20. 

n 

N =-1025.7 

1 

-1025.7 

2 

Sandwich 

Compr. 

. 30 

.025 

20000 

20. 

X 

Bl 

. 30 

.025 

50000 

.25 

60. 

3 

N =-1025.7 

n 

-1025.7 

Sandwich 

Compr. 

n 

. 30 

.025 

20000 

20. 

■1 

Si 

. 30 

.25 

B 

B 

20. 

1 

N =-11975 

16 

-11975 

Isotropic 

Compr. 

ffij 

. 30 

.25 

B 

B 

20. 

JJS 

. 30 

.25 

B 

m 

20. 

wm 

N =27931 

16 

27963 

Isotropic 

Shear 

. 30 

.25 

1 

B 

20. 

mi 

. 30 

1/6 

B 

1 

17. 

wm 

N =11455 

16 

11468 

Isotropic 

Shear 

. 30 

1/6 

B 

B 

17. 

Hi 

xy 

.25 

.025 

4359. 

.25 

10. 

10 

N =1344.3 

B 

1349.2 

Sandwich 

Shear 

. 25 

.025 

4359. 

40. 

xy 

v(% 

.25 

5496. 

. 1875 

19 

10 

N =1252.7 

B 

1260.5 

Sandwich 

Shear 

.25 

BS 

5496. 

xy 

m 

lO.fcxlO6 

.25 

.005 

4359. 

.25 

24. 

11 

NxB*  =258.7 

10 

256. 1 

Sandwich 

Bending 

10.6x10 

.25 

.005 

4359. 

40. 

■ 

lO.ZxlO6 

. 30 

.005 

.225 

48. 

D 

N _ =119.  3 

1 

119.8 

10 

Sandwich 

Bending 

10.2x10 

. 30 

.005 

48. 

IBr  Hill 

. 30 

.0625 

B 

40. 

B 

N -,=117. 1 

16 

116.7 

)1 

Isotropic 

Bending 

. 30 

.0625 

D 

M 

60. 

M 

Maximum  value  along  edge  is  given  here;  i.  e.  , N^  = Nxg(l-2y/b). 


25 


Comparison  of  the  calculations  for  pure  edgewise  shear  loads  show 
much  larger  discrepancies  than  for  bending  or  compression.  This  is  to 
be  expected,  since  convergence  to  the  critical  shear  load  is  quite  slow  due 
to  the  marked  difference  between  the  forms  of  the  assumed  displacements 
and  the  true  mode  shapes.  Slightly  lower  values  are  obtained  for  each 
example  by  the  present  method.  These  results  appear  to  be  more  accurate 
than  the  previously  published  values,  since  more  terms  in  the  displacement 
series  (Equations  3.  1)  are  considered,  and  because  the  solution  does  not 
depend  upon  discarding  various  terms  from  the  complete  series,  as  has 
been  done  in  Reference  7. 

4.  2 TYPICAL  RESULTS  FOR  COMBINED  LOADINGS 

Critical  loads  calculations  for  typical  sandwich  panels  subjected  to 
combined  loadings  are  presented  in  Figures  4 and  5.  The  physical  data  for 
each  of  the  three  panels  considered  is  listed  in  Table  2.  Panel  1 is  square, 
with  balanced  facings  and  equal  shear  moduli  in  the  core.  Panel  2 is  similar 
to  Panel  1,  but  rectangular  (a/b  = .556).  The  third  panel  studied  is  rectangular 
with  unequal  core  shear  moduli. 

The  effect  of  edgewise  shear  loading  upon  the  magnitude  of  the  critical 

compressive  loads  is  plotted  in  Figure  4.  Here  N = N for  all  points 

xo  yo 

calculated,  while  the  magnitude  of  the  applied  shear  load  is  increased  rela- 
tive to  the  compressive  forces.  Results  are  given  in  terms  of  the  ratio 

N 

_xc r _ Actual  critical  compressive  load 

N * Critical  compressive  load  for  N =0 
xc  r r xy 

and  the  relative  magnitude  of  the  shear  and  compressive  stress,  N /N  . 

xy  x 

Figure  5 depicts  the  effect  of  the  addition  of  edgewise  bending  moments 

to  an  applied  biaxial  loading.  In  each  case  the  compressive  loadings  N , 

xo 

N are  equal.  A shear  load  is  also  applied  such  that  N =2N  . Various 

yo  xy  xo 


fixer* 

1.00 

0.90 

0.80 

0.70 

0.60 

0.50 

0.40 

0.30 

0.20 

0.10 


-3-2-101234 

^xB^xo 

Figure  5.  Effect  of  an  Edgewise  Bending  Load  Upon  Critical 
Compressive  and  Shear  Stresses. 
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TABLE  2:  PHYSICAL  DATA  FOR  COMBINED -LOADS  EXAMPLES 


♦ 


Panel  1 

Panel  2 

Panel  3 

C=  E 

E1 

9. 5xl06 

10. 6xl06 

10. 2x1 0^ 

E2 

9. 5xl06 

10. 6xl06 

10. 2x1 0^ 

S 

. 021 

. 050 

. 030 

lz 

.021 

. 050 

. 030 

v. 

. 30 

. 30 

. 28 

V 

2 

. 30 

. 30 

. 28 

G 

19000. 

30000. 

8000. 

cxz 

G 

19000. 

30000. 

22000. 

cyz 

t 

. 181 

. 25 

. 370 

c 

a 

20. 

20. 

32. 

b 

20. 

36. 

64. 
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proportions  of  edgewise  bending  stresses  are  then  added  to  this  state  of 
stress,  and  the  resulting  loads  observed.  Results  are  expressed  in  terms 
of  the  ratio 


r 


I 


xoc  r Actual  critical  compressive  load 

N * Critical  compressive  load  for  N _=0 
xoc  r xB 


and  the  relative  magnitude  of  the  bending  and  direct  compressive  stresses, 


/N 


xo* 


4.  3 COMPARISON  WITH  INTERACTION  FORMULAS 

In  order  to  facilitate  the  analysis  of  stability  for  cases  of  combined 
loadings,  various  interaction  formulas  have  been  proposed.  A number  of 
examples  are  given  in  Reference  1.  The  purpose  of  any  interaction  formula 
is  to  account  for  the  combined  effects  of  various  applied  loads  by  com- 
parison of  the  intensity  of  each  load  to  the  critical  value  of  the  same  load 
acting  alone.  Each  such  ratio  is  weighted  by  means  of  an  exponent, 
associated  with  the  particular  type  of  loading;  for  example, 

el  e2 

Rj  + R^  = const.  , 

where 

R.  = N./N. 
i i icr 

The  e^  are  constant  exponents,  and  the  index  i is  associated  with  the  type 
of  loading. 

The  use  of  interaction  formulas  is  attractive  since  any  number  of 
load  cases  may  be  analyzed  once  the  individual  critical  loads  are  deter- 
mined. However,  such  formulas  obviously  represent  approximations  to  an 
exact  analysis  of  stability  under  combined  loadings.  Interaction  formulas 
for  the  cases  of  edgewise  compression  combined  with  edgewise  shear, 
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edgewise  compression  combined  with  edgewise  bending,  and  edgewise  shear 
combined  with  edgewise  bending  have  been  evaluated  to  a limited  extent 
during  the  course  of  this  study. 


4.  3.  1 Interaction  Formula  for  Combined  Edgewise  Shear  and 
Compression 

The  interaction  formula  for  axial  compression  and  edgewise 

shear  is^ 

N /N  + (N  /N  )2  = 1.  (4.  1) 

x xcr  xy  xycr 

When  biaxial  compression  is  present.  Equation  4.  1 is  generalized  to  the 

, 8 

form 

N /N  + N /N  + (N  /N  )2  = 1.  (4.2) 

x xcr  y ycr  xy  xycr 

The  left-hand  side  of  Equation  4.  2 has  been  evaluated  (see 
Table  3)  for  Panels  1 and  2 of  Table  2,  for  various  combinations  of  biaxial 
compression  and  shear  loads.  The  formula  is  reasonably  accurate  for  most 
of  the  cases  considered.  For  the  square  panel  (Panel  1)  the  maximum  error 
observed  is  .28%.  For  the  rectangular  panel  (Panel  2)  this  error  increases 
to  more  than  11%.  It  is  important  to  realize  that  a value  of  the  left-hand 
side  of  Equation  4.2  which  is  greater  than  one  represents  a conservative 
estimate  of  the  critical  loads.  For  all  cases  included  in  Table  3,  the  inter- 
action formula  is  conservative,  with  the  exception  of  one  example  (Case  2-7) 
involving  shear  combined  with  both  tension  and  compression  loads. 

Equation  4.  1 has  also  been  examined  for  various  proportions 
of  compression  and  shear  loads.  The  error  in  Equation  4.  1 is  plotted  in 
Figure  6 as  a function  of  these  proportions  (a  positive  error  corresponds 
to  a conservative  result,  as  noted  previously).  Results  for  the  square  panel 
were  again  very  good.  Cases  tested  for  the  rectangular  panel  showed  a 
maximum  error  of  approximately  4.  5%.  Figure  7 shows  the  compression - 
shear  interaction  curves  predicted  by  the  present  method  in  relation  to  a 
plot  of  the  interaction  formula  (Equation  4.  1). 
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TABLE  3:  COMPARISON  WITH  INTERACTION  FORMULA 

FOR  COMBINED  BIAXIAL  COMPRESSION  AND 
EDGEWISE  SHEAR 


n 

xo 

n 

xy 

n 

yo 

N 

xcr 

N 

xycr 

N 

yc  r 

Equation 
4.  2 

i 

i 

I 

-421.40 

±881. 33 

-421. 40 

i 

i 

1 

-200. 09 

-200. 09 

-200.  09 

1. 0019 

i 

2 

1 

-177. 16 

-354. 32 

-177. 16 

1. 0024 

2 

2 

1 

-257.  44 

-257.44 

-128. 72 

1. 0017 

1 

3 

2 

-118. 11 

-354. 33 

-236. 22 

1. 0025 

3 

1 

0 

-411.45 

-137. 15 

0. 

1. 0006 

1 

5 

2 

- 97.66 

-488. 30 

-195.  32 

1. 0022 

2 

6 

3 

-133. 98 

-401.94 

-200. 97 

1.0028 

1 

1 

1 

1 

2 

1 

2 

1 

5 

1 

-1037. 24 

±3453.94 

-2331. 08 

-759. 70 

-759. 70 

-759 • 70 

1. 1125 

-683.49 

-1366.98 

-683.  49 

1. 0942 

-851.82 

-851.82 

-425.91 

1. 0648 

-696.  10 

-1740.25 

-348. 05 

1. 0743 

-693. 27 

-1617.63 

-462.  18 

1.0860 

-553.76 

-2353. 48 

-138.44 

1.0576 

-921.8 

I 

1 

1843.6 

460.  90 

.9759 

00  1 


4.  3.  2 Interaction  Formula  for  Combined  Edgewise  Bending  and 
Compression 


The  interaction  formula  for  combined  edgewise  bending  and 
axial  compression  loads  is  given  by* 


N /N 
xo  xocr 


(N  /N  )3/2  = 1 . 
xB  xBcr 


(4.3) 


Here  N ^ is  a pure-bending  load  as  shown  in  Figure  1.  Calculated  values 
of  the  left-hand  side  of  the  formula  are  listed  in  Table  4,  and  the  corre- 
sponding errors  are  plotted  in  Figure  8 as  a function  of  the  relative  pro- 
portions of  the  two  loads.  Results  for  the  rectangular  panel  using 
Equation  4.  3 are  reasonably  good.  The  maximum  error  (5.  3%)  occurs  for 
a moderately  large  bending  load.  For  the  square  panel,  Equation  4.3  gives 
accurate  results  only  when  the  bending  load  is  quite  small  compared  to  the 
compressive  force.  When  the  bending  load  predominates,  the  error  be- 
comes as  large  as  18%.  Again,  the  interaction  formula  provides  a con- 
servative estimate  of  the  critical  load  in  each  case  considered.  The  same 
results  are  compared  graphically  in  Figure  9 with  the  shape  of  the  inter- 
action curve  (Equation  4.  3). 


4.  3.  3 Interaction  Formula  for  Combined  Edgewise  Bending  and 
Shear 

The  interaction  formula  proposed*  for  the  evaluation  of  cri- 
tical combined  bending  and  shear  loads  is 


(N  / N ) 
xB  xBcr 


+ (N  /N 

xy  xycr 


(4.4) 


Values  of  the  left-hand  side  of  Equation  4.  4,  as  given  in  Table  5,  indicate 
nonconservative  buckling  predictions  by  the  interaction  formula.  Errors 
for  the  square  panel  considered  are  less  than  8%.  Those  for  the  rectangular 
panel  are  as  high  as  13.  5%.  These  values  are  plotted  in  Figure  10  for 
various  relative  proportions  of  bending  and  shear.  Critical  load  combina- 
tions are  compared  to  the  shape  of  the  interaction  curve  in  Figure  11, 
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TABLE  4:  COMPARISON  WITH  INTERACTION  FORMULA  FOR 

COMBINED  EDGEWISE  BENDING  AND  AXLAL 
COMPRESSION 


Case 

Panel 

n 

xo 

11  14 

xB 

N 

xocr 

— 

N n 
xBc  r 

Equation 
4.  3 

1 

i 

i 

———————— 

-421.40 

±2387. 55 

1-1 

i 

±i 

-410.67 

410.  67 

1. 0459 

1-2 

2 

±i 

-418. 60 

209. 30 

1.0193 

1-3 

3 

±i 

-420. 14 

140  °5 

1.0112 

1-4 

6 

±i 

-421. 08 

70.  18 

1. 0043 

1-5 

1 

±2 

-384. 17 

768. 34 

1.0942 

1-6 

1 

±4 

-320. 17 

1280. 69 

1. 1526 

1-7 

1 

±5 

-291. 51 

1457. 56 

1.  1688 

1-8 

1 

±6 

-266. 43 

1598. 58 

1.  1801 

1 

1 

-1037. 24 

±4436.90 

2-1 

2 

1 

±1 

-972. 09 

972. 09 

1.0  397 

2-2 

2 

±1 

-1019. 04 

509. 52 

1.0214 

2-3 

3 

±1 

-1028.97 

342. 99 

1. 0135 

2-4 

6 

±1 

-1035. 14 

172. 52 

1. 0056 

2-5 

1 

±2 

- 847. 36 

1694. 71 

1. 0530 

2-6 

1 

±4 

- 636.  13 

2544. 50 

1.0476 

2-7 

1 

± 5 

- 560. 71 

2803.  54 

1.0429 

2-8 

1 

±6 

- 500. 21 

3001. 28 

1. 0386 

in 


o 


o' 

m 


dOddB 


w 

V 

U 
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COMPRESSION1 


TABLE  5:  COMPARISON  WITH  INTERACTION  FORMULA  FOR 
COMBINED  EDGEWISE  BENDING  AND  SHEAR 


Case 

Panel 

"xB 

— 

n 

xy 

N » 
xB 

N 

xy 

Equation 
4.  4 

1 

i 

±2387. 55 

±883.  78 

1-1 

1 

1 

i 

815. 29 

815.29 

.9676 

1-2 

- 1 

i 

815. 29 

815. 29 

.9676 

1-3 

2 

i 

1370.  34 

685.  17 

.9305 

1-4 

4 

i 

1900. 11 

475.03 

. 9223 

1-5 

1 

3 

291.69 

875. 07 

.9953 

1-6 

1 

5 

176. 12 

880. 61 

.9983 

1-7 

1 

8 

110. 32 

882. 53 

.9993 

1-8 

1 

10 

88.  30 

882. 98 

.9996 

1 

1 

±4436. 90 

±3453. 94 

2-1 

2 

1 

1 

2534. 14 

2534. 14 

. 8645 

2-2 

- 1 

1 

-2534.  14 

2534. 14 

. 8645 

2-3 

2 

1 

3534. 07 

1767.  04 

. 8962 

2-4 

4 

1 

4131. 24 

1032. 81 

.9564 

2-5 

1 

3 

1074. 01 

3222. 02 

.9288 

2-6 

1 

5 

668. 70 

3343.49 

. 9598 

2-7 

1 

8 

425. 15 

3401.19 

. 9789 

2-8 

1 

10 

341. 79 

3417.89 

.9852 

39 


BENDING  : 


40 


Combined  Edgewise  Shear  and  Bending. 


4.4  DISCUSSION 


Various  calculations  concerning  the  stability  of  a sandwich  panel  under 
combined  modes  of  loading  and  the  effectiveness  of  typical  interaction 
formulas  have  been  presented  in  this  section.  These  results  indicate  the 
character  of  some  effects  due  to  combined  edge  loadings,  and  illustrate  the 
versatility  of  the  analysis  procedure  described  herein. 

The  data  given  in  Paragraph  4.  2 indicates  that  small  amounts  of 
bending  or  shear  loading  added  to  large  compressive  forces  are  usually  not 
critically  important  (see  Figures  4 and  5).  This  can  be  deduced  also  from 
the  large  magnitudes  of  the  critical  bending  and  shear  stresses  as  compared 
with  those  for  compression  (Tables  3 and  4).  The  effect  of  a bending  load 
upon  critical  biaxial  stresses  shows  a rather  wide  variation  for  different 
panel  geometry  and  stiffness,  in  contrast  to  the  influence  of  a shear  load 
upon  critical  compressive  stresses. 

Results  given  concerning  the  various  interaction  formulas  indicate 
that  the  formulas  proposed*  for  combined  edgewise  shear  and  biaxial 
compression  (Equation  4.2)  and  for  combined  compression  and  bending 
(Equation  4.  3)  may  be  of  some  use  for  design  purposes,  since  both  are 
conservative  for  all  cases  considered  here.  For  the  case  of  combined  ten- 
sion, compression  and  shear,  an  example  has  been  considered  (Table  3, 

Case  2-7)  for  which  the  interaction  formula  (Equation  4.2)  yields  a non- 
conservative result.  The  shear -uniaxial  compression  formula  (Equation  4.1) 
generally  gives  good  results  and  is  conservative.  The  equation  for  bending 
and  compression  should  be  used  with  some  caution  in  the  event  that  bending 
loads  become  large.  The  formula  for  combined  edgewise  bending  and  shear 
(Equation  4.4)  appears  to  be  unsuitable  for  the  determination  of  buckling 
loads,  since  the  error  incurred  by  its  use  is  not  only  large,  but  noncon- 
servative. 
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SECTION  5 

I 

SUMMARY  AND  CONCLUSIONS 

An  analysis  procedure  has  been  described  for  the  prediction  of  the 
general  instability  of  flat,  simply- supported , rectangular  sandwich  panels 
loaded  by  arbitrary  combinations  of  biaxial  edgewise  compression,  biaxial 
edgewise  bending,  and  edgewise  shear  forces.  The  mathematical  model  is 
formulated  in  terms  of  the  total  potential  energy,  which  is  composed  of  the 
internal  elastic  strain  energy  of  the  sandwich  plate  and  the  potentials  of  the 
applied  loadings.  The  Ritz  method  is  used  to  transform  the  energy  functional 
in  terms  of  continuous  field  variables  into  a quadratic  function  of  certain 
discrete  parameters.  Application  of  the  principle  of  minimum  potential 
energy  results  in  a generalized  discrete  eigenvalue  problem,  which  can  be 
solved  by  standard  methods.  The  finite  number  of  eigenvalues  obtained  in 
the  solution  approximate  the  lowest  of  the  infinite  number  of  critical  loads, 
while  the  eigenvectors  approximate  the  corresponding  buckled  mode  shapes. 

The  accuracy  of  the  analysis  has  been  demonstrated  by  comparison 
with  results  available  in  the  literature  for  a number  of  examples.  A limited 
study  has  been  made  of  the  applicability  of  typical  interaction  formulas  which 
are  commonly  used  in  sandwich  design.  Comparisons  of  the  interaction 
formulas  with  numerical  results  generated  from  the  present  analysis  indi- 
cates that  certain  of  these  formulas  can  yield  substantially  nonconservative 
results  for  some  combinations  of  geometry  and  loading,  and  should  be  used 
with  caution. 

A computer  program  has  been  developed  for  the  implementation  of  the 
combined- loads  instability  analysis.  The  program  is  operable  both  inter- 
actively and  by  batch  mode  processing.  A detailed  description  of  the  com- 

, . , . . 17 

puter  program  and  its  usage  are  contained  in  a companion  report. 
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The  analysis  presented  has  been  demonstrated  to  be  a useful  tool  for 
the  stability  analysis  of  a practical  class  of  sandwich  panels.  However,  a 
number  of  extensions  could  be  undertaken  to  apply  the  analysis  to  more 
general  situations.  Remaining  within  the  confines  of  the  linear  formulation, 
extensions  of  the  present  method  to  include  layered  composite  faces,  elastic 
edge  members,  and  more  general  panel  shapes  should  be  given  consideration. 
The  analysis  can  be  generalized  to  the  cases  of  multicore  sandwich  and  panels 
having  normal  deformations,  at  some  expense  in  the  number  of  degrees  of 
freedom  to  be  considered.  The  consideration  of  boundary  conditions  other 
than  simple  supports  are  accomodated  by  altering  the  assumed-mode  func- 
tions to  be  used  in  the  Ritz  approximations.  Hov.  .ver,  the  computational 
effort  is  increased  for  more  general  support  conditions,  since  the  normal 
displacement  does  not  uncouple  from  the  other  displacement  variables,  as  is 
the  case  for  simply  - supported  boundaries  (see  Equations  3.6  and  3.7). 

In  conclusion,  the  sandwich  instability  analysis  described  in  this  report 
represents  a useful  and  accurate  tool  for  the  design  and  analysis  of  light- 
weight, high-performance  structural  components.  In  addition,  the  method 
reported  provides  a suitable  starting  point  for  the  consideration  of  many  more 
general  types  of  sandwich  construction. 
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